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A Monte Carlo generator has been constructed to simulate the reaction e + e~ — > 7 + 2n and 7 + 4tt where the 
photon is assumed to be observed in the detector. Predictions are presented for cms energies of lGeV, 3GeV 
and lOGeV, corresponding to the energies of DA$NE, BEBC and of B-meson factories. The event rates are 
sufficiently high to allow for a precise measurement of R(Q 2 ) in the region of Q between approximately lGeV and 
2.5GeV. Estimates for the kinematic breaking of isospin relations between different channels as a consequence of 
the 7T - 7T° mass difference are given. 



1. Introduction 

The precise determination of the cross section 
for electron positron annihilation into hadrons 
over a large energy range is one of the important 
tasks of current particle physics. The results are 
relevant for the analysis of electroweak precision 
measurements which are affected by the running 
of the electromagnetic coupling from the Thomp- 
son limit up to Mz- Also the prediction for the 
anomalous magnetic moment of the muon de- 
pends critically on these data. Last not least the 
measurement of the energy dependence of R(s) is 
one of the gold plated tests of QCD and allows 
for a precise determination of the strong coupling 
constant. 

Depending on the energy region different tech- 
niques for the measurement of R(s) have been 
applied up to date. At low energies, say from the 
two pion threshold up to roughly two GeV, exclu- 
sive channels are collected separately. For higher 
energies inclusive measurements start to become 
dominant. For energies below m r isospin invari- 
ance and CVC have traditionally been used to 
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predict r decays from electron positron annihila- 
tion 0-§. Clearly this strate gy can be inverted 
|^] , pending irreducible uncertainties from isospin 
violation and radiative corrections 

To cover a large range of energies, results from 
many different experiments and colliders have to 
be combined, and energy scans have to be per- 
formed to obtain the full energy dependence. An 
attractive alternative is provided by the upcom- 
ing and B— meson factories which operate at 
large luminosities, albeit at fixed energies. Events 
with radiated tagged photons give access to a 
measurement of R over the full range of energies, 
from threshold up to the CMS energy of the col- 
lider. For events with tagged photons the invari- 
ant mass of the recoiling hadronic system is fixed 
by the photon energy which provides an impor- 
tant kinematic constraint. 

To arrive at reliable predictions including angu- 
lar and energy cutoffs as employed by realistic ex- 
periments, a Monte Carlo generator is indispens- 
able. For hadronic states with invariant masses 
below two or even three GeV it is desirable to sim- 
ulate the individual exclusive channels with two, 
three up to six mesons, i.e. pions, kaons, etas etc. 
which requires a fairly detailed parametrization 
of the various form factors. 
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In principle initial and final state radiation 
would be required for the complete simulation. 
Such a program has been constructed for the two 
pion case . There it is demonstrated that suit- 
ably chosen configurations, namely those with 
hard photons at small angles relative to the beam 
and well separated from the pions, are dominated 
by initial state radiation. In fact, this separation 
is possible |J even when operating the <fi factory 
DA$NE on top of the <j> resonance where direct 
radiative <fi decays cannot be ignored. 

In a recent paper || this project has been con- 
tinued with the construction of a generator for the 
radiative production of the four pion final state, 
including the lu(— > 37r)7r channel. This mode con- 
tributes a large fraction of the rate with invariant 
masses between one and two GeV. The energy 
region between 1.5 GeV and 2.5 GeV is difficult 
to access directly with current electron positron 
colliders. At the same time the experimental un- 
certainties are relatively large. This motivates the 
special effort devoted to this range. 

The Monte Carlo program is constructed in a 
modular form such that the parametrization of 
the hadronic matrix element can easily be re- 
placed by a more elaborate version. Different final 
states with three, four or five pions or kaons can 
be included. The present parametrization of the 
hadronic matrix element follows closely the form 
suggested in ]Tc| ] , correcting only some minor de- 
ficiencies. The four pion amplitude is assumed to 
be dominated by pi — > ixa\ plus a direct coupling 
pi — > pirir and exhibits the proper behavior in the 
chiral limit. 

2. The radiative return 

Hard photons observed at small angles relative 
to the electron or positron beam and at the same 
time well separated from charged particles in the 
final state can be used to reduce the effective 
center of mass energy at electron positron col- 
liders. As shown in Table |l[ the event rates are 
fairly high, of O(10 6 ) both at DAPHNE and at 
B-meson factories. Performing a detailed analy- 
sis of the angular and energy distributions for the 
77r + 7r~ final state it has been shown that initial 
and final state radiation can be reasonably well 



separated @||. For the four pion case we there- 
fore restrict the discussion to initial state radia- 
tion only. The amplitude for the production of 
an arbitrary hadronic final state 
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involves the matrix element of the hadronic cur- 
rent 

jr = »•,«») 

ee <%!),.., h(q n )\Jr(0)\0> (2) 

which has to be parameterized by form factors 
to be discussed below. For the two pion case the 
amplitude is determined by only one function, the 
pion form factor F%^. 

JT'** = {qt - Qu) F 2 , (0 2 ) (3) 

The matrix element for the four pion case will be 
discussed below. After integrating the hadronic 
tensor H av over the hadronic phase space one gets 



(4) 
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{Q»Qu - g^Q 2 ) R(Q 2 



Gtt 

where R(Q 2 ) is a(e + e~ —> hadrons)/a po i n t- The 
additional integration of the differential cross sec- 
tion over the photon angles (the azimuthal angle 
is integrated over the full range and the polar an- 
gle within Q m in < < ir — 9 m in ) leads to the 
differential distribution 
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2 ) 



cos 9„ 



(5) 



which can be used to calculate the event rate 
observed for realistic photon energy and angular 
cuts (see Tab.l). 

3. Isospin relations 

The emphasis of || was towards final states 
consisting of four pions and a photon. Be- 
fore entering a discussion of a model depen- 
dent parametrization of the form factors the 
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Event rates 


Collider 




Integrated luminosity, fb 1 


"mm — ° 


"min — ' 


#min — 10° 


DA$NE 


1.02 


1 


13 • 10 B 


12 • 10 B 


10 • 10 y 


B — factory 


10.6 


100 


4- 10 B 


3.5 • 10 b 


3 ■ 10 b 


B — factory 


10.6 


100 


2.7- 10 a 


2.3 • 10 B 


2.0 • 10 6 



Table 1 

Estimated number of radiative events e + e~ — * hadrons + 7 for different center of mass energies. In 
the first two rows hadrons stands for ir + ir~ and the minimal photon energy is 0.1 GeV. The third row 
is obtained for a continuum contribution in the region 2 GeV < \J Q 2 < 3.7 GeV assuming a constant 
R = 2.4. 



constraints from isospin invariance must be re- 
called. They relate the amplitudes of the e + e _ — > 
2ir + 2ir~ and e + e _ — > ir + ir~2ir° processes and 
those for r decays into 7r _ 37r° and tt + 2tt^tt q . The 
amplitude for r decay 

G 

M T = -j= cos 6 C v (p v ) 7 Q (1 - 75) u (p T ) J~ (6) 



leads to the differential distribution 
dT 

dO 2 ~ 



(7) 



2 r e ^(i-^V(i + 2 ^W 



with 



J ~-J v * d$ n (Q; qi ,...,q n ) 



(8) 



3?r 



(Q»Qu-g»uQ 2 ) R T (Q 2 ) 



Note the relative factor of 2 between the defini- 
tions in Eq.(§) and Eq.(g). 

Final states with an even number of pions are 
produced through the isospin one part of the elec- 
tromagnetic current only, whence 



V2 J« m (27r) = J-(27r) (9) 
and R(Q 2 ) = R T (Q 2 ) for two pion final states. 



A similar relation for the four pion final state 
is easily obtained H||: 

(7r+7r"<7r 2 |J^|0> = M P i,p 2 ,P + ,p~) 
<7r+7r+ 7^1^10) = 

Jn(pi>P2,Pl>Pl) + J »(Pl:P2:Pt:Pl) 
+ J /*(P2 ) Pi: Pi >P2 ) + J iiipl:Pl:Pi,P2~) 

(7r-7r°^7r§|J-|0) = J M (p 2 ,p 3 ,P" ,Pi)+ 

Jn(Pl,P3,P~,P2) + Jn( P l,P2,P~,P3) 

(7rr7rj7r + 7r°|J-|0) = 

J»{p + :P2:Pt:P°) + Mp + , Pi, P2, P°) , (10) 

which connects r decay and electron positron an- 
nihilation (J em = -Lj 3 ). The (in Q 2 ) differ- 
ential rates are then consistent with the familiar 
relations between t decays and e + e~ annihilation 
into four pions: 

/T(-OOO) = l R (+ + — ) (11) 
R T (--+Q) = ^R{+ + — )+i?(+-oo) 



4. The 7T ± -7r° mass difference 

All the relations obtained in the previous sec- 
tion are strictly applicable only in case all pions in 
the final states have the same mass, which is obvi- 
ously not true. The relatively large (about 3.6%) 
tt^-it mass difference will affect the R{Q 2 ) 
R T (Q 2 ) relations even if Eq.@ and Eq.@ still 
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only P A modified 


(3 6 and T p modified 


0.28 < E [GeV] < 0.81 
Aa(M z ) 


1.0163 
1.0116 


1.0088 
1.0027 


0.32 < E [GeV] < 1.777 
a M 
Aa(Mz) 


1.0130 
1.0096 


1.0058 
1.0016 



Table 2 

Kinematic correction factors for the predictions of and Aa from r data. 



hold. Let us retain these relations between the 
amplitudes and incorporate the tt^-t: mass dif- 
ference in the phase space considerations. This 
should give at least an indication of the size of 
these "kinematic" isospin violations. 

To estimate these effects for the two-pion case 
assume that the cross section a(e + e~ — > ir + ir~) 
is well measured and the squared form factor 



\F(Q 2 )\ 2 = 



3Q 2 

ita 2 (3 3 



a(e+e- ^tt+tt-)(Q 2 ) , (12) 



(j3„ = (l - 4to 2 ._/<3 2 ) 1 ^ then used to predict 

the — > i/ r 7r _ 7r° decay rate The size of the cor- 
rections depends critically on the details of the 
assumptions. 

If the form factor and the form of the current 

J e ; n (2n) = V2 J-(2n) = (q hfi - q 2 ^) F(Q 2 ) (13) 

remain unchanged, the integral for the r rate re- 
ceives an S-wave contribution which can be elim- 
inated by replacing Eq.fO) by 



V2J-(2n) 



~ 92,/. 



(14) 



Q 2 



Q.(qi - q 2 ) F(Q 2 ) 



Numerically the contribution to the integral of 
this latter term is tiny - nevertheless we shall 
adopt the second choice. 

The introduction of the tt ± — 7r° mass differ- 
ence raises the prediction for the r decay rate 
by 0.86%. It seems, however, plausible, that the 
energy dependent width of the p- meson, which 



is present in the form factor, has to be modi- 
fied accordingly, leading to an effective increase 
of T p by 0.74%. The two effects nearly com- 
pensate in the integral. Hence the relation be- 
tween t~ — > v T ir~ir° partial decay width and the 
e + e~ — * 7r + 7r~ cross section would only be cor- 
rected by 0.06%. However, a sizable Q 2 depen- 
dence of the ratio of the two spectral functions 
(V(Q 2 ) ~ \F{Q 2 )\ 2 (3 3 ) is expected, with a re- 
duction approximately 0.74 % close to the peak 
of the p resonance and enhancements at the tails. 

At present, however, r data provide an impor- 
tant input for the prediction of the QED coupling 
at the scale of Mz and the hadronic contribution 
to g-2 |j 11 1. Let us, for the moment, assume 
that the aforementioned kinematic effects are in- 
deed present. If the e + e~ cross section is deduced 
from t data through Eq. ( ^6|) 



1 dT(T- 



ZA-7T 



'7T°) _ 3COS 2 C 



dQ 2 

2X 2 



2na 2 m 2 . 



1-^) [1 + 2^1 *(e 



mz 



(15) 



7T + 7T ) 



and the contributions to g-2 and a(Mz) are evalu- 
ated without kinematic corrections of phase space 
and form factor the former are overestimated by 
0.58% and 0.16% respectively (Tab.2). 

The situation is more complicated for the four 
pion case. To estimate the size of the effect the 

quantities R( hO) etc. were evaluated once 

assuming that all masses are equal to to_ and 
once taking the real masses. The corrections to 
the integrated quantities modify the integrated 
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version of the Eq. ■ 

— !— IY-OOO) 
1.050 V ' 

r(--+o) 

1.024 V ; 



accordingly: 
1 



r(+ + 
r(+ + 



(16) 
(17) 



1 



+ L046 r ( + -°°) 



5. The hadronic current 

Employing isospin relations it is thus sufficient 
to construct the hadronic current for the (H — 00) 
mode only. The basic building block of the cur- 
rent contains a part built on the assumption of 
ai vector dominance plus an lo exchange contri- 
bution. By adding an /o contribution one indeed 
recovers the proper chiral limit p2| . The three 
contributions are depicted schematically in Fig.[l] 
and described in detail in 

To test the physical predictions of the cur- 
rent, let us start with the r decay branching ra- 
tios, summarized in Tab. 3. The agreement be- 
tween the predictions of the current || and the 
data is satisfactory for the t~ — > j/ r 7r _ 37r° de- 
cay mode. Comparing however the results for the 
t~ — > v T 27r~ tt + n° and t~ — > h> T iT~L) (7r~7r + 7r°) 
modes it seems that the uj part of the current 
does not represent the data well, even if the total 
branching ratio for the r~ — > v t 2it~ it + it® decay 
mode agrees with the data. 

In |^] also the predictions for the e + e~ — > 
2ir + 2ir~ and e + e~ — » tt + tt^2tt q cross sections 
are compared with data, and satisfactory agree- 
ment between the Monte Carlo and the data is 
observed. 

6. The Monte Carlo program 

The idea behind the structure of the Monte 
Carlo program is to allow for a simple addition of 
new final state modes into the program and for a 
simple replacement of the current (s) of the exist- 
ing modes. The program thus exhibits a modu- 
lar structure. For the generation of the four mo- 
menta of the mesons no sophisticated method of a 
variance reduction was applied. This slows down 
the generation, but could be accounted for if a 
faster Monte Carlo generator would be required. 



It has, however, the advantage of being universal 
and no change of the variance reduction method 
is required with each modification of the hadronic 
current. The process to be simulated by the pro- 
gram in its final stage is e + e _ — > 7 + hadrons 
with an exclusive description of final states, even 
if till now only ir + ir~, 2tt°tt + tt~ and 2tt + 2tt^ 
hadronic final states are implemented. The LL ra- 
diative QED corrections were taken into account 
using structure function method as developed in 
|lg| ] and limited to the initial emission only. In 
fact the program can run in one of two modes 
(chosen by a user) one with collinear radiation 
and one without it. Hard large angle photon emis- 
sion is limited to initial state radiation, which is 
justified by where it was demonstrated for the 
7r + 7r~ hadronic state that the contribution from 
the final state emission as well as the initial-final 
state interference can be reduced to a negligible 
level by applying suitable cuts. 

Let us now discuss the cuts, which reduce the 
contribution of final state radiation to a negligible 
level. We recall that in @ it was shown that the 
following set of angular cuts 



cutsl : (18) 
(7° < 7 < 20° or 160° < 7 < 173°) 
and 30° <6 n < 150° , 



( 



'7 



is the photon (pion) polar angle) ful- 



fils this requirement for the e H 



7 cross 



section at DAPHNE. It reduces, however, the ob- 
served cross section significantly. This starts to 
become dramatic, when one runs at energies well 
above 1 GeV. The following set of cuts 



cuts2 : 



(19) 



(7° < 6> 7 < 20° and 30° < 9 n < 173°) 
or (160° < 7 < 173° and 7° < 6^ < 150°) 

also reduces the contribution from final state ra- 
diation to a negligible level due to the fact that 
the pions and photon are well separated as in the 
previous case. 

At the same time the cross section reduction is 
much smaller, especially for higher beam energies 



G 




Figure 1. Diagrams contributing to the hadronic current. 



Mode 


|1C] 


present model 


experiment 


Br (r~ -> ^ T 27r"7r+7r ) 


3.11% 


4.33% 


4.20(8)% 


Br (t~ — ► z/ T 7r _ ti; (7r _ 7r + 7r )) 


1.20% 


1.48% 


1.73(11)% 


Br (r - -> i/ T 7r-37T u ) 


0.98% 


1.14% 


1.08(10)% 



Table 3 

Branching ratios of r decay modes. Results of flo|| and the present current are compared to experimental 
data. 
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% [nb/GeV]_ 



1 1 

no cuts v 
with cutsl ❖ 
with cuts2 t 




1.5 2 
/(?[GeV] 



2.5 3 



Figure 2. The differential e + e~ — > 2ir + 2ir~^ 
cross sections for a beam energy of 1.5 GeV and 
a minimal photon energy of 0.1 GeV for different 
cuts as described in the text. 



and higher energies of the observed photons. The 
effect of the two sets of cuts on the cross sections 
with four charged pions in the final state is pre- 
sented in Fig|| for 2E beam = 3 GeV. "No cuts" 
corresponds to 7° < Q 1 < 173° and not cuts on 
the pions. For higher beam energies the effect of 
the cross section reduction is much stronger and 
at 2 Ei, eam = 10 GeV the cross sections for the 
cuts specified in Eq.(pj|) is reduced almost to zero. 
For the cuts specified in Eq. ( p0| ) the reduction re- 
mains tolerable (Table §). 

One may even restrict photon and pions detec- 
tion angles to the central region, e.g. 25° < 6* 7 < 
155° and 30° < 9 n < 150° respectively if a min- 
imal angle of 20° between photon and charged 
and neutral pions is required in order to suppress 
final state radiation and to clearly separate neu- 
tral pions and the photon. With this cut one 
obtains (y/s = 10 GeV and C = 100 fb" 1 ) a rate 
of 1.17 • 10 5 events with 2tt + 2tt--i and 1.53 • 10 5 
events with 27r°7r + 7r~7. 

Additional collinear emission always present in 



the real experiment reduces slightly the cross sec- 
tions. Its actual size depends on the cuts on the 
invariant mass of the 4-7T + 7 system. The ef- 
fect is similar for different energies and for both 
charge modes. 

7. Summary 

As an alternative to a direct measurement of 
the cross section for e + e~ — > hadrons at the 
relevant energy one may use initial state radi- 
ation to reduce the effective energy of electron 
positron colliders, exploiting the large luminosity 
of "factories" and accessing thus a continuum of 
hadronic final states. 

With this motivation a Monte Carlo genera- 
tor has been constructed to simulate the reaction 
e + e~ — » 7 + 27r and 7 + 47r, where the photon 
is assumed to be observed in the detector. Once 
more accurate data become available, the modu- 
lar structure of the program will allow for modi- 
fication or replacement of the hadronic current in 
a simple way. Additional collinear photon radia- 
tion has been incorporated with the technique of 
structure functions. 

Predictions are presented for cms energies of 
lGeV, 3GeV and lOGeV, corresponding to the 
energies of DA$NE, BEBC and of £?-meson fac- 
tories. Even after applying realistic cuts the event 
rates are sufficiently high to allow for a precise 
measurement of R{Q 2 ) in the region of Q between 
approximately lGeV and 2.5GeV. 

Acknowledgments. The author would like to 
thank H. Czyz for the pleasant collaboration on 
the topics discussed in this talk. Work supported 
by BMBF under grant BMBF-057KA92P. 

REFERENCES 

1. H.B. Thacker and J.J. Sakurai, Phys. Lett. 
36B (1971) 103. 

2. P. Tsai, Phys. Rev. D4 (1971), 2821. 

3. F.J. Gilman and D.H. Miller, Phys. Rev. D17 
(1978) 1846; F.J. Gilman and S.H. Rhie Phys. 
Rev. D31 (1985) 1066. 

4. J.H. Kiihn, A. Santamaria, Zeit. f. Phys. C48 
(1990) 445. 



8 







Event rates 
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10 GeV 
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